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Abstract: Lehmann's ideas on concepts of dependence have had a profound 
effect on mathematical theory of reliability. The aim of this paper is two-fold. 
The first is to show how the notion of a "hazard potential" can provide an 
explanation for the cause of dependence between life-times. The second is to 
propose a general framework under which two currently discussed issues in re- 
liability and in survival analysis involving interdependent stochastic processes, 
can be meaningfully addressed via the notion of a hazard potential. The first 
issue pertains to the failure of an item in a dynamic setting under multiple 
interdependent risks. The second pertains to assessing an item's life length 
in the presence of observable surrogates or markers. Here again the setting is 
dynamic and the role of the marker is akin to that of a leading indicator in 
multiple time series. 



1. Preamble: Impact of Lehmann's work on reliability 

Erich Lehmann's work on non-parametrics has had a conceptual impact on reliabil- 
ity and life-testing. Here two commonly encountered themes, one of which bears his 
name, encapsulate the essence of the impact. These are: the notion of a Lehmann 
Alternative, and his exposition on Concepts of Dependence. The former (see 
Lehmann [4]) comes into play in the context of accelerated life testing, wherein a 
Lehmann alternative is essentially a model for accelerating failure. The latter (see 
Lehmann [5]) has spawned a large body of literature pertaining to the reliability 
of complex systems with interdependent component lifetimes. Lehmann's original 
ideas on characterizing the nature of dependence has helped us better articulate 
the effect of failures that are causal or cascading, and the consequences of lifetimes 
that exhibit a negative correlation. The aim of this paper is to propose a frame- 
work that has been inspired by (though not directly related to) Lehmann's work on 
dependence. The point of view that we adopt here is "dynamic" , in the sense that 
what is of relevance are dependent stochastic processes. We focus on two scenarios, 
one pertaining to competing risks, a topic of interest in survival analysis, and the 
other pertaining to degradation and its markers, a topic of interest to those work- 
ing in reliability. To set the stage for our development we start with an overview of 
the notion of a hazard potential, an entity which helps us better conceptualize the 
process of failure and the cause of interdependent lifetimes. 
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2. Introduction: The hazard potential 

Let T denote the time to failure of a unit that is scheduled to operate in some 
specified static environment. Let h(t) be the hazard rate function of the survival 
function of T, namely, P(T > t),t > 0. Let H(t) = J Q h(u)du, be the cumulative 
hazard function at t;H(t) is increasing in t. With h(t),t > specified, it is well 
known that 

Pr(T > t; h(t),t > 0) = exp(-H(t)). 

Consider now an exponentially distributed random variable X, with scale parameter 
A, A > 0. Then for some H(t) > 0, 

Pr(X > H(t)\X = 1) = exp(-fl"(t)); 

thus 

(2.1) Pr(T > t; h(t), i > 0) = exp(-H(t)) = Pr{X > H(t)\X = 1). 

The right hand side of the above equation says that the item in question will 
fail when its cumulative hazard H(t) crosses a threshold X, where X has a unit 
exponential distribution. Singpurwalla [11] calls X the Hazard Potential of the 
item, and interprets it as an unknown resource that the item is endowed with at 
inception. Furthermore, H(t) is interpreted as the amount of resource consumed 
at time t, and h{t) is the rate at which that resource gets consumed. Looking at 
the failure process in terms of an endowed and a consumed resource enables us to 
characterize an environment as being normal when H(t) = t, and as being acceler- 
ated (decelerated) when H(t) > (<) t. More importantly, with X interpreted as an 
unknown resource, we are able to interpret dependent lifetimes as the consequence 
of dependent hazard potentials, the later being a manifestation of commonalities 
of design, manufacture, or genetic make-up. Thus one way to generate dependent 
lifetimes, say T\ and T 2 is to start with a bivariate distribution (Xi,X 2 ) whose 
marginal distributions are exponential with scale parameter one, and which is not 
the product of exponential marginals. The details are in Singpurwalla [11]. 

When the environment is dynamic, the rate at which an item's resource gets 
consumed is random. Thus h(t);t > is better described as a stochastic process, 
and consequently, so is H(t),t > 0. Since H(t) is increasing in t, the cumulative 
hazard process {H(t);t > 0} is a continuous increasing process, and the item 
fails when this process hits a random threshold X, the item's hazard potential. 
Candidate stochastic processes for {H(t); t > 0} are proposed in the reference given 
above, and the nature of the resulting lifetimes described therein. Noteworthy are 
an increasing Levy process, and the maxima of a Wiener process. 

In what follows we show how the notion of a hazard potential serves as a unifying 
platform for describing the competing risk phenomenon and the phenomenon of 
failure due to ageing or degradation in the presence of a marker (or a bio marker) 
such as crack size (or a CD4 cell count). 

3. Dependent competing risks and competing risk processes 

By "competing risks" one generally means failure due to agents that presumably 
compete with each other for an item's lifetime. The traditional model that has 
been used for describing the competing risk phenomenon has been the reliability of 
a series system whose component lifetimes are independent or dependent. The idea 
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here is that since the failure of any component of the system leads to the failure 
of the system, the system experiences multiple risks, each risk leading to failure. 
Thus if Tj denotes the lifetime of component i, i = 1, . . . , k, say, then the cause 
of system failure is that component whose lifetime is smallest of the k lifetimes. 
Consequently, if T denotes a system's lifetime, then 

(3.1) Pr(T > t) = P(JTi(t) <X U ..., H k (t) < X k ), 

where Xi is the hazard potential of the i-th component, and Hi(t) its cumula- 
tive hazard (or the risk to component i) at time t. If the X^'s are assumed to be 
independent (a simplifying assumption), then (3.1) leads to the result that 

(3.2) Pr(T > t) = eM~(Hi(t) + ■■■ + H k (t))], 

suggesting an additivity of cumulative hazard functions, or equivalently, an ad- 
ditivity of the risks. Were the X^s assumed dependent, then the nature of their 
dependence will dictate the manner in which the risks combine. Thus for example 
if for some 6, < 6 < 1, we suppose that 

Pr(-Xi > Xi,X 2 > x 2 \6) = cxp(— x\ — x 2 — 6x^X2), 

namely one of Gumbcl's bivariate exponential distributions, then 

Pr(T > t\6) = exp[-(ffi(t) + H 2 (t) + OH^H^t))]. 

The cumulative hazards (or equivalently, the risks) are no longer additive. 

The series system model discussed above has also been used to describe the 
failure of a single item that experiences several failure causing agents that compete 
with each other. However, we question this line of reasoning because a single item 
posseses only one unknown resource. Thus the Xi, . . . ,X k of the series system model 
should be replaced by a single X, where X\ = X 2 = ■ ■ ■ = X k = X (in probability). 
To set the stage for the single item case, suppose that the item experiences k 
agents, say C\, . . . ,C k , where an agent is seen as a cause of failure; for example, 
the consumption of fatty foods. Let Hi(t) be the consequence of agent d, were d 
be the only agent acting on the item. Then under the simultaneous action by all of 
the k agents the item's survival function 

PT(T>t;h 1 (t),...,h k (t)) 

(3.3) =P(H 1 (t)<X,...,H k (t)<X) 

= exp(-max(# 1 (i),...,# fe (*))). 

Here again, the cumulative hazards are not additive. 

Taking a clue from the fact that dependent hazard potentials lead us to a 

non- additivity of the cumulative hazard functions, we observe that the condition 

p p p p p 
Xi = X 2 = ■ ■ ■ = X k = X (where X\ = X 2 denotes that Xiand X 2 are equal in 

probability) implies that X\ , . . . , X k are totally positively dependent, in the sense 

of Lehmann (1966). Thus (3.2) and (3.3) can be combined to claim that in general, 

under the series system model for competing risks, P(T > t) can be bounded as 

k 

(3.4) exp{-J2 H i(t))< p ( T >t) <exp(-max(JTi(t),...,iTk(i))). 

1 

Whereas (3.4) above may be known, our argument leading up to it could be new. 
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3.1. Competing risk processes 

The prevailing view of what constitutes dependent competing risks entails a con- 
sideration of dependent component lifetimes in the series system model mentioned 
above. By contrast, our position on a proper framework for describing dependent 
competing risks is different. Since it is the Hi(t)'s that encapsulate the notion 
of risk, dependent competing risks should entail interdependence between Hi(t)'s, 
i = 1, . . . , k. This would require that the £/j(£)'s be random, and a way to do so 
is to assume that each {Hi(t);t > 0} is a stochastic process; we call this a com- 
peting risk process. The item fails when any one of the {Hi(t);t > 0} processes 
first hits the items hazard potential X. To incorporate interdependence between 
the Hi(t)'s, we conceptualize a fe-variate process {Hi(t), . . . , Hk(t); t > 0}, that we 
call a dependent competing risk process. Since i?i(i)'s are increasing in i, one 
possible choice for each {Hi(t); t > 0} could be a Brownian Maximum Process. That 
is Hi(t) = sup 0<s<t {Wi(s); s > 0}, where {Wi(s);s > 0} is a standard Brownian 
motion process. Dependence between the iSi(i)'s can be induced via a dependence 
between the {Wi(s);s > 0} processes. Thus for example, in the bivariate case, if 
p denotes the correlation between two standard Brownian motion processes, then 



Pr(T > t) 



P(Hi(t) < x,H 2 (t) < x)e- x dx 



and it can be shown (details omitted) that, 



(3.5) 



Pr(T > t) 



So So CX P 


(a 2 +b 2 -2pab) 
2*(1-^) 


dadb 


Jo Jo ex P 


[u 2 -\-v 2 ~ 2puv} 
2*(1-P 2 ) 


dudv 



Another possibility again for the case of k = 2, is to assume that {Hi(t); t > 0} 
is some non- negative, non-decreasing, right-continuous process, but that {H 2 (t)-, 
t > 0} has a sample path which is an impulse function of the form H2(t) = for all 
t < t*, and that H 2 (t*) = oo for some t* > 0, where the rate of occurrence of the 
impulse at time t depends on H\(t). The process {H2(t);t > 0} can be identified 
with some sort of a traumatic event that competes with the process {H\(t);t > 0} 
for the lifetime of the item. In the absence of trauma the item fails when the 
process {Hi(t);t > 0} hits the item's hazard potential. This scenario parallels the 
one considered by Lemoine and Wenocur [6], albeit in a context that is different 
from ours. By assuming that the probability of occurrence of an impulse in the time 
interval [t,t + h), given that H\(t) = ui, is 1 — cxp(— uih), Lemoine and Wenocur 
[6] have shown that for X = x, the probability of survival of an item to time t is of 
the form: 



(3-6) 



Pr(T >t) = E 



exp( / Hi(s)dfl) J [0iX) (ffi(t)) 



where Ia{*) is the indicator of a set A, and the expectation is with respect to the 
distribution of the process {H\{t);t > 0}. As a special case, when {H\(t);t > 0} is 
a gamma process (see Singpurwalla [10]), and x is infinite, so that /[o )00 ) (Hi(t)) = 1 
for Hi(t) > 0, the above equation takes the form 



(3.7) 



Pr(T > t) = exp(-(l + t) log(l + 1) + t). 
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The closed form result of (3.7) suffers from the disadvantage of having the effect of 
the hazard potential de facto nullified. The more realistic case of (3.6) will call for 
numerical or simulation based approaches. These remain to be done; our aim here 
has been to give some flavor of the possibilities. 

4. Biomarkers and degradation processes 

A topic of current interest in both reliability and survival analysis pertains to assess- 
ing lifetimes based on observable surrogates, such as crack length, and biomarkers 
like CD4 cell counts. Here again the hazard potential provides a unified perspective 
for looking at the interplay between the unobservable failure causing phenomenon, 
and an observable surrogate. It is an assumed dependence between the above two 
processes that makes this interplay possible. 

To engineers (cf. Bogdanoff and Kozin [1]) degradation is the irreversible accu- 
mulation of damage throughout life that leads to failure. The term "damage" is 
not defined; however it is claimed that damage manifests itself via surrogates such 
as cracks, corrosion, measured wear, etc. Similarly, in the biosciences, the notion 
of "ageing" pertains to a unit's position in a state space wherein the probabilities 
of failure are greater than in a former position. Ageing manifests itself in terms 
of biomedical and physical difficulties experienced by individuals and other such 
biomarkers. 

With the above as background, our proposal here is to conceptualize ageing and 
degradation as unobservable constructs (or latent variables) that serve to describe 
a process that results in failure. These constructs can be seen as the cause of ob- 
servable surrogates like cracks, corrosion, and biomarkers such as CD4 cell counts. 
This modelling viewpoint is not in keeping with the work on degradation modelling 
by Doksum [3] and the several references therein. The prevailing view is that degra- 
dation is an observable phenomenon that reveals itself in the guise of crack length 
and CD4 cell counts. The item fails when the observable phenomenon hits some 
threshold whose nature is not specified. Whereas this may be meaningful in some 
cases, a more general view is to separate the observable and the unobservable and 
to attribute failure as a consequence of the behavior of the unobservable. 

To mathematically describe the cause and effect phenomenon of degradation (or 
ageing) and the observables that it spawns, we view the (unobservable) cumula- 
tive hazard function as degradation, or ageing, and the biomarker as an observ- 
able process that is influenced by the former. The item fails when the cumulative 
hazard function hits the item's hazard potential X, where X has exponential (1) 
distribution. With the above in mind we introduce the degradation process as 
a bivariate stochastic process {H(t),Z(t),t > 0}, with H{t) representing the un- 
observable degradation, and Z(t) an observable marker. Whereas H(t) is required 
to be non-decreasing, there is no such requirement on Z(t). For the marker to be 
useful as a predictor of failure, it is necessary that H(t) and Z(t) be related to each 
other. One way to achieve this linkage is via a Markov Additive Process (cf. Cinlar 
[2]) wherein {Z(t);t > 0} is a Markov process and {H(t);t > 0} is an increasing 
Levy process whose parameters depend on the state of the {Z(t);t > 0} process. 
The ramifications of this set-up need to be explored. 

Another possibility, and one that we are able to develop here in some detail (see 
Section 5), is to describe {Z(t); t > 0} by a Wiener process (cracks do heal and CD4 
cell counts do fluctuate), and the unobservable degradation process {H(t);t > 0} 
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by a Wiener Maximum Process, namely, 

(4.1) H(t) = sup {Z(s);s> 0}. 

0<s<t 

What makes the topic of analyzing degradation processes attractive is not just 
the modeling part; the statistical and computational issues that the set-up creates 
are quite challenging. Since {Z(t); t > 0} is an observable process, how may one use 
observations on this process until some time, say t*, to make inferences about the 
process of interest H(t), for t > t*l In other words, how does one assess Pr(T > 
t\{Z(s); < s < t* < t}) , where T is an item's time to failure? Furthermore, as is 
often the case, the process {Z(s);s > 0} cannot be monitored continuously. Rather, 
what one is able to do is observe {Z(s);s > 0} at k discrete time points and use 
these as a basis for inference about Pr(T > t\{Z(s);0 < s < t* < t}). These and 
other matters are discussed next in Section 5, which could be viewed as a prototype 
of what else is possible using other models for degradation. 

5. Inference under a Wiener maximum process for degradation 

We start with some preliminaries about a Wiener process and its hitting time to a 
threshold. The notation used here is adopted from Doksum [3]. 

5. 1 . Hitting time of a Wiener maximum process to a random threshold 

Let Z t denote an observable marker process {Z(t);t > 0}, and H t an unobservable 
degradation process {H(t);t > 0}. The relationship between these two processes 
is prescribed by (4.1). Suppose that Z t is described by a Wiener process with a 
drift parameter r\ and a diffusion parameter a 2 > 0. That is, Z(0) = and Z t has 
independent increments. Also, for any t > 0, Z(t) has a Gaussian distribution with 
E{Z{t)) = rft, and for any < h < t 2 , Var[Z(t 2 ) - Z(ti_)} = (t 2 - h)a 2 . Let T x 
denote the first time at which Z t crosses a threshold x > 0; that is, T x is the hitting 
time of Z t to x. Then, when rj = 0, 

Pr (Z{t) > x) = Pr(Z(i) > x\T x < t) Pr(T x < t) 

+ Pr (Z(t) > x\T x > t) Pr (T x > t) , 

(5.2) Pr (T x <t) = 2Pr(Z(t)>x). 

This is because Pr(Z(t) > x\T x < t) can be set to 1/2, and the second term on 
the right hand side of (5.1) is zero. When Z(t) has a Gaussian distribution with 
mean rjt and variance a 2 t, Pr(Z(f) > x) can be similarly obtained, and thence 

def 

Pr(T x < t) = F x (t\n,cr). Specifically it can be seen that 

where [i = x/n and A = x 2 /a 2 . The distribution F x is the Inverse Gaussian 
Distribution (/G-Distribution) with parameters a and A, where a = E(T X ) and 
\[i 2 =Var(T :E ). Observe that when i] = 0, both E{T X ) and V&r(T x ) are infinite, and 
thus for any meaningful description of a marker process via a Wiener process, the 
drift parameter i] needs to be greater than zero. 
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The probability density of F x at t takes the form: 

A (t-rf 



(5.4) /,■(.'!'/■") = \/ exp 



2u 2 



for t,u,X> 0. 

We now turn attention to H tl the process of interest. We first note that because 
of (4.1), H(0) = 0, and H(t) is non-decreasing in t; this is what was required of 
H t . An item experiencing the process H t fails when H t first crosses a threshold 
X, where X is unknown. However, our uncertainty about X is described by an 
exponential distribution with probability density f(x) = e~ x . Let T denote the 
time to failure of the item in question. Then, following the line of reasoning leading 
to (5.1), we would have, in the case of n = 0, 

Pr(T <t) = 2Pr(H(t) > x). 

Furthermore, because of (4.1), the hitting time of H t to a random threshold X will 
coincide with T x , the hitting time of Z t (with 77 > 0) to X. Consequently, 

Pr (T <t) = Pr(T x <t) = Pr(T x < t\X = x)f(x)dx 

Jo 

Pr{T x <t)e- x dx= F x (t\r),a) e - X dx. 



roc 



Rewriting F x (t\i],a) in terms of the marker process parameters 77 and a, and treating 
these parameters as known, we have 

Pv(T<t\ V ,a) = f F(t\ V ,a) 

<"» -fK^-^) + *(-^-^) 

x exp ^x ^— I — 1^ dx, 

as our assessment of an item's time to failure with 77 and er assumed known. It is 
convenient to summarize the above development as follows 

Theorem 5.1. The time to failure T of an item experiencing failure due to ageing 
or degradation described by a Wiener Maximum Process with a drift parameter 
n > 0, and a diffusion parameter a 2 > 0, has the distribution function F(t\n,a) 
which is a location mixture of Inverse Gaussian Distributions. This distribution 
function, which is also the hitting time of the process to an exponential (1) random 
threshold, is given by (5.5). 

In Figure 1 we illustrate the behavior of the /G-Distribution function F x (t), 
for x = 1,2,3,4, and 5, when r\ = <r = 1, and superimpose on these a plot of 
F(t\n = a = 1) to show the effect of averaging the threshold x. As can be expected, 
averaging makes the 5-shapedness of the distribution functions less pronounced. 



5.2. Assessing lifetimes using surrogate (biomarker) data 



The material leading up to Theorem 5.1 is based on the thesis that r\ and a 2 are 
known. In actuality, they are of course unknown. Thus, besides the hazard potential 
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Fig 1. The IG- Distribution with thresholds x = 1, . . . , 5 and the averaged IG- Distribution. 

X, the rj and a 2 constitute the unknowns in our set-up. To assess rj and a 2 we may 
use prior information, and when available, data on the underlying processes Z t and 
H t . The prior on X is an exponential distribution with scale one, and this prior 
can also be updated using process data. In the remainder of this section, we focus 
attention on the case of a single item and describe the nature of the data that can 
be collected on it. We then outline an overall plan for incorporating these data into 
our analyses. 

In Section 5.3 we give details about the inferential steps. The scenario of observ- 
ing several items to failure in order to predict the lifetime of a future item will not 
be discussed. 

In principle, we have assumed that H t is an unobservable process. This is cer- 
tainly true in our particular case when the observable marker process Z t cannot be 
continuously monitored. Thus it is not possible to collect data on H t . Contrast our 
scenario to that of Doksum [3], Lu and Meeker [7], and Lu, Meeker and Escobar 
[8], who assume that degradation is an observable process and who use data on 
degradation to predict an item's lifetime. We assume that it is the surrogate (or 
the biomarker) process Z t that is observable, but only prior to T, the item's failure 
time. In some cases we may be able to observe Z t at t = T, but doing so in the case 
of a single item would be futile, since our aim is to assess an unobserved T. Data 
on Z t will certainly provide information about r\ and cr 2 , but also about X; this is 
because for any t < T, we know that X > Z(t). Thus, as claimed by Nair [9], data 
on (the observable surrogates of) degradation helps sharpen lifetime assessments, 
because a knowledge of 77, a 2 and X translates to a knowledge of T. 

It is often the case - at least we assume so - that Z t cannot be continuously 
monitored, so that observations on Z t could be had only at times < t\ < t^ < ■ ■ ■ < 
tk < T, yielding Z = (Z(ti), . . . , Z{tk)) as data. Furthermore, based on Z(tk) 7 we 
are able to assert that X > Zitk). This means that our updated uncertainty about 
X will be encapsulated by a shifted exponential distribution with scale parameter 
one, and a location (or shift) parameter Z(tk). 

Thus for an item experiencing failure due to degradation, whose marker process 
yields Z as data, our aim will be to assess the item's residual life (T — tk). That is, 
for any u > 0, we need to know Pr(T > tk + u; Z) = Pr(T > tk + u; T > tk), and 
this under a certain assumption (cf. Singpurwalla [12]) is tantamount to knowing 




(5.6) 



Pr(T > t k + it) 
Pr(T > t k ) 



On Competing risk and degradation processes 



237 



for < u < oo. To assess the two quantities in the above ratio, we need to consider 
the quantity Pr(T > t; Z), for some t > 0. Let n(r], <r 2 ,x; Z) encapsulate our uncer- 
tainty about 77, a 2 and X in the light of the data Z. In Section 5.3 we describe our 
approach for assessing ^(77, cr 2 , x; Z). Now 

(5.7) Pr(T>t;Z)=/ Pr(T > % a 2 , x; Z)n(r], a 2 , x; Z)(dr 1 )(d<T 2 ){dx) 

J r),a 2 ,x 

= I Pv{T x >t\r],<T 2 )ir(ri,a 2 ,x;Z)(dri){d<T 2 )(dx) 

(5.8) = f F x {t\r],a)^a 2 ,x;Z){dvi){da 2 ){dx), 

where F x (t\r],a) is the /G-Distribution of (5.3). 

Implicit to going from (5.7) to (5.8) is the assumption that the event (T > t) 
is independent of Z given 77, a 2 and X. In Section 5.3 we will propose that 77 be 
allowed to vary between a and b; also, a 2 > 0, and having observed Z(tk), it is clear 
that x must be greater than Z(tk). Consequently, (5.8) gets written as 

f-b f-OC r-OO 

(5.9) Pr(T>t;Z)=/ / / F x (t\ V , <t)k( V , a 2 , x; Z)(d V )(da 2 )(dx), 

Ja JO JZ(t k ) 

and the above can be used to obtain Pr(T > rj fe + u; Z) and Pr(T > t^; Z). Once 
these are obtained, we are able to assess the residual life Pr(T > + u\T > tk), 
for u > 0. 

We now turn our attention to describing a Bayesian approach specifying 7r(r;, a 2 , 
x;Z). 

5.3. Assessing the posterior distribution of 77, cr 2 and X 

The purpose of this section is to describe an approach for assessing 7r(?7, cr 2 , x; Z), 
the posterior distribution of the unknowns in our set-up. For this, we start by 
supposing that Z is an unknown and consider the quantity ir(i],a 2 ,x\ Z). This is 
done to legitimize the ensuing simplifications. By the multiplication rule, and using 
obvious notation 

Trfa, cr 2 , x\Z) = 7T\ (77, a 2 \X, Z)n 2 (X\Z). 
It makes sense to suppose that 77 and a 2 do not depend on X; thus 

(5.10) n(r], n 2 ,x\Z) = n^n, <j 2 \Z)tt 2 (X\Z). 
However, Z is an observed quantity. Thus (5.10) needs to be recast as: 

(5.11) 7r(77,cr 2 ,a;;Z) =7ri(77,cr 2 ;Z)7r 2 (X;Z). 

Regarding the quantity TT2(X;Z), the only information that Z provides about 
X is that X > Z(tk). Thus 7r 2 (X;Z) becomes ^{X; Z(tk))- We may now invoke 
Bayes' law on n2(X; Z(tk)) and using the facts that the prior on X is an exponential 
(1) distribution on (0, 00), obtain the result that the posterior of X is also an 
exponential (1) distribution, but on (Z(tk), 00). That is, tt2(X; Z(tk)) is a shifted 
exponential distribution of the form exp(— (x — Z(tk))), for x > Z(tk)- 

Turning attention to the quantity 7Ti(?7, cr 2 ; Z) we note, invoking Bayes' law, that 



(5.12) 



71-1(77, cr 2 ; Z) oc C(n, cr 2 ; Z)tt*(77, cr 2 ), 
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where C(r], cr 2 ; Z) is the likelihood of 77 and a 2 with Z fixed, and tt* (77, cr 2 ) our prior 
on 77 and a 2 . In what follows we discuss the nature of the likelihood and the prior. 

The Likelihood of 77 and a 2 

Let Y x = Z(h), Y 2 = (Z(t 2 ) - Z(h)), . . . , Y k = (Z(t k ) - Z(t k -i)), and 8l = 
£1 , S2 = ti — ti, . . . , Sfe = tk — tk—i- Because the Wiener process has independent 



increments, the y,'s are independent. Also, j/j ~ N(r]Si,a Sj), 



1, . . . , k, where 



N(fi,l; 2 ) denotes a Gaussian distribution with mean \x and variance £ 2 . Thus, the 
joint density of the t/j's, i = 1, . . . , k, which is useful for writing out a likelihood of 
7] and a 2 , will be of the form 



m - ys t 

(T 2 Si 



where </> denotes a standard Gaussian probability density function. As a consequence 
of the above, the likelihood of r/ and a 2 with y = (j/i, . . . , y^) fixed, can be written 
as: 



(5.13) 



>C(r/, ( r 2 ;y) = n 



exp 



Vi - V s i 



The Prior on rj and cr 2 

Turning attention to 7r* (77, cr 2 ), the prior on 77 and cr 2 , it seems reasonable to suppose 
that 77 and cr 2 are not independent. It makes sense to suppose that the fluctuations 
of Zt depend on the trend 77. The larger the 77, the bigger the cr 2 , so long as there is 
a constraint on the value of 77. If 77 is not constrained the marker will take negative 
values. Thus, we need to consider, in obvious notation 

(5.14) 7T*(77,Cr 2 ) - 7r*(cr 2 |77)7T*(7i). 

Since 77 can take values in (0, 00), and since 77 = tan 9 - see Figure 2-9 must 
take values in (0, tt/2). 

To impose a constraint on 77, we may suppose that 9 has a translated beta density 
on (a, b), where < a < b < ir/2. That is, 9 = a + (b — a)W, where W has a beta 
distribution on (0, 1). For example, a could be ir/8 and b could be 3tt/8. Note that 
were 9 assumed to be uniform over (0, n/2), then 77 will have a density of the form 
2/[-7r(l + i] 2 )] - which is a folded Cauchy. 




Fig 2. Relationship between Zt and rj. 
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The choice of ir*(a 2 \r]) is trickier. The usual approach in such situations is to 
opt for natural conjugacy. Accordingly, we suppose that ip ^— has the prior 



where v is a parameter of the prior. 

Note that E(tp\-q,v) = 77/(1/ — 2), and so ip — ct 2 increases with 77, and r\ is 
constrained over a and b. Thus a constraint on a 2 as well. 

To pin down the parameter v, we anchor on time t = 1, and note that since 
E(Z\) = r\ and Var(Zi) = a 2 = ip, a should be such that Act should not exceed 
77 for some A = 1,2,3,...; otherwise Z\ will become negative. With A = 3, 77 = 
3ct and so ip = a 2 = r] 2 /9. Thus v should be such that E(cr 2 \r],v) w ?7 2 /9. But 
E{a 2 \rj 1 v) = rj/(v — 2), and therefore by setting r\j(y — 2) = ry 2 /9, we would have 

v = 9/r/ + 2. In general, were we to set 77 = Act, v = A 2 /i] + 2, for A = 1, 2, 

Consequently, u/2 + l = (A 2 /t? + 2)/2 + 1 = A 2 /2r7 + 2, and thus 



would be our prior of ct 2 , conditioned on ip, and A = 1,2, . .., serving as a prior 
parameter. Values of A can be used to explore sensitivity to the prior. 

This completes our discussion on choosing priors for the parameters of a Wiener 
process model for Z t . All the necessary ingredients for implementing (5.9) are now 
at hand. This will have to be done numerically; it does not appear to pose major 
obstacles. We are currently working on this matter using both simulated and real 
data. 

6. Conclusion 

Our aim here was to describe how Lehmann's original ideas on (positive) depen- 
dence framed in the context of non-parametrics have been germane to reliability 
and survival analysis, and even so in the context of survival dynamics. The notion 
of a hazard potential has been the "hook" via which we can attribute the cause 
of dependence, and also to develop a framework for an appreciation of competing 
risks and degradation. The hazard potential provides a platform through which the 
above can be discussed in a unified manner. Our platform pertains to the hitting 
times of stochastic processes to a random threshold. With degradation modeling, 
the unobservable cumulative hazard function is seen as the metric of degradation 
(as opposed to an observable, like crack growth) and when modeling competing 
risks, the cumulative hazard is interpreted as a risk. Our goal here was not to solve 
any definitive problem with real data; rather, it was to propose a way of looking at 
two commonly encountered problems in reliability and survival analysis, problems 
that have been well discussed, but which have not as yet been recognized as having 
a common framework. The material of Section 5 is purely illustrative; it shows what 
is possible when one has access to real data. We are currently persuing the details 
underlying the several avenues and possibilities that have been outlined here. 
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